London theory across superconducting phase transitions: apphcation to UPts 
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For multi-component superconductors, it is known tliat the presence of symmetry breaking fields 
can lead to multiple superconducting phase transitions. This phenomenon is best illustrated in 
UPts. Motivated by recent small angle neutron scattering experiments on the vortex state of UPts, 
the London theory in the vicinity of such phase transitions is determined. It is found that the form 
of this London theory is in general quite different than that for conventional superconductors. This 
is due to the existence of a diverging correlation length associated with these phase transitions. One 
striking consequence is that non-trivial vortex lattices exist arbitrarily close to Ifci . It is also found 
that the penetration depth develops a novel temperature dependence and anisotropy. Results of this 
theory for UPta are derived. Possible applications to Celus, Ui_a;Tha;Bei3, electron doped cuprate 
superconductors, Sr2Ru04, and MgCNia are also discussed. 

One of the qualitatively new phenomena associated with unconventional superconductors is the existence of multiple 
degrees of freedom in the Cooper pair wavefunction This leads naturally to the possibility of phase transitions 

between superconducting states of different symmetry. A well known example is UPta , which exhibits a phase diagram 
with three superconducting phases, each described by an order parameter with different symmetry |^-|^. Recently, 
using small angle neutron scattering (SANS), Huxley et al. observed a realignment of the flux-line lattice as the 
phase boundary between two different phases of UPts is crossed This measurement provides the motivation 
to derive the London theory in the vicinity of such transitions (see Ref. Q] for a recent study of this problem). 
Conceptually this raises the question: does the phase transition lead to non-trivial consequences for the low field 
vortex phases? Here we address this question and find that the London theory indeed develops a novel form near 
the transition. The resulting theory profoundly changes the vortex phase diagram from that expected from the usual 
London theory. One result is the appearance of a rectangular flux lattice structure at the lower critical field (-ffci)- 
Such a lattice structure is not thought to be possible in the dilute vortex limit near Hd ■ Physically, this arises because 
the coherence length associated with the second transition diverges and thus becomes an important length scale near 
Hci- The diverging coherence length allows a non-trivial anisotropy to exist arbitrarily close to Hd- This anisotropy 
dictates the appearance of the resulting rectangular vortex phase. We show that such behavior should exist in UPta. 
Possible appHcations to Ui_a;Tha;Bei3, Celns, electron doped cuprate superconductors, MgCNia, and Sr2Ru04 are 
also discussed. 
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FIG. 1. Superconducting phase diagram of UPta 
for fields near Hd- 



Temperature 

We are interested in the vortex phases near the ^ to _B phase transition 



Ginzhurg Landau Free energy 

To demonstrate the physics described above the superconducting state of UPta will be considered. The supercon- 
ducting phase diagram of UPts exhibits three distinct superconducting phases (see Figure 1). There have been many 
theoretical interpretations of this phase diagram (see Ref. |^ for an overview). However, it appears that only those 
theories in which a doubly degenerate order parameter has its degeneracy broken by a symmetry breaking field are 
consistent with the recent SANS measurements of Huxley et al. In particular, it was observed that the low-field 
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flux-lattice in the ^-phase is not aligned with any of the crystallographic axes when the magnetic field is applied 
along the six-fold symmetry axis. |^. This observation implies the existence of a two-component order parameter 
since it cannot occur for any of the single-component representations of Dqh ■ The two most popular variants of 
this class of models are based on the i?2« and the Eig representations of the hexagonal point group ||^,^ . The origin 
of the symmetry breaking field is taken to be antiferromagnetic order |^,^. However, the results are valid for any 
weak symmetry breaking field that reduces the Dqh symmetry to D2h- We consider the £'211 and the Eig theories 
with an applied magnetic field along the c-axis and ignore spatial variations along this axis. The resulting free energy 
for the i?2u or Eig representations of Dqh is 

/ = - a\rj\^ + 7(77+77* + 77_77* ) + \[3M^ + [32\v+?\v-? 
K|D?7+|2 + |Dr;_|2 

+ k[{D^r^+y{D+ri^) + c.c.) (1) 

where a = ao{Tc — T), Dj = Vj — ^Aj, h = V x A, and A is the vector potential, and D± = Dj.± iDy. This free 
energy neglects a term hz{\r]^\'^ — |?7_P) which will not change any of the main physical results. This term has been 
argued to be small Furthermore, there is a coupling of the symmetry breaking field to the gradient terms in the free 
energy which has not been included Q . This coupling is required to give rise to the observed phase diagram at high 
fields 1^, but does not result in any significant changes in the low-field regime discussed in this paper. Experiments 
constrain the phenomenological parameters for UPta. By fitting to the size of the specific heat jumps the value of 
P2/ Pi = 0.75. From the splitting of the two transition temperatures |^ 7/010 = 17 mK (this gives Tc^a = 500 mK and 
Tc^B — 440 mK). These parameters will be used throughout this paper. The crucial parameter that drives the main 
results of this work is k. The magnitude of k/k, has been the subject of debate |^,||. Sauls has argued that k/n, ||] 
is small while Joynt ||^ has argued that k/n may be as large as 1. In the following we will perform a perturbation in 
k/ K, and show that even a moderate k/n has profound consequences on the vortex phase diagram for T « Tc,b- 

The goal of this work is to develop the London theory for the field along the six-fold symmetry axis as the AtoB 
phase transition is crossed. This problem has been addressed recently by Champel and Mineev in Ref. 0. The result 
of this study is that the vortex lattice will be a stretched hexagonal lattice near Hd- This was found by determining 
the linear combination of 77+ and ?7_ that stabilizes the homogeneous ground state (this linear combination will be 
denoted below by 7/"+) and by examining the spatial variations of this order parameter. Such an analysis will not 
result in the physics discussed in the introduction. This physics is found by including the order parameter component 
orthogonal to that which minimizes the homogeneous minimizing free energy (denoted by ip-). The component ■(/)_ 
is induced by spatial variations of "0+ through the k term in the free energy. As will be shown below, this coupling 
qualitatively changes physics of the London theory from that found by setting ~ 0. 

The homogeneous free energy stabilizes [rj+^rj^) = (1, 1) near Tc+ (the A phase) and at a lower temperature there 
is a second transition to the B phase, a state {r]+,ri-) = |r7| [cos((/9/2), sin((^/2)] where if ^ 7r/2 and varies with 
temperature. This transition occurs when (52 > 0. This free energy will be used to determine the London theory with 
the lowest order non-local correction (similar to what was done for the superconducting state of Sr2Ru04 ||^). To 
implement this the (?7+,77_) is rotated through an angle ip/2 to give a basis {tp+Tip-)- The angle Lp is chosen so that 
(-0+, ^Z"-) oc (1, 0) is the stable solution for all temperatures. For 7 > 0, this requires that = 7r/2 in the A phase and 
that siiiip — and = a in the B phase. The rotated free energy is 

/ = - a\iP\^+j[cosip{ij+ij*- + V'-V4) + sin^'dV^+l' - + ^AlV^I" + 

i/32[(|r/+P + 177-12)2 - [cos ^(V'+^l + sin ^(IV^+P ~ IV- P)]'] 
A^dD^+P + IDV^.p) 

+ '^k[[{D^i;+r{D+^^) - (i?_7/^_)*(Z?+V+) + c.c] + cos(^[(i?_V^+)*(Z?+V_) + {D^^^r{D+4>+) + c.c] 

+ sin(^[(i?_V^_)*(i?+V^_) - {D^ib+nD+ib+) + c.c]] (2) 

where cc denotes complex conjugate. If k = then the two components ip^ and ip- are decoupled, implying that the 
solution in the vortex phase near Hd is described by tp- = 0. Once k becomes non-zero, this term induces a non-zero 
in the vortex phase. We treat the term in k up to second order perturbatively. All terms in 0{^p^), 0{kip'^) and 
higher are ignored in the free energy which results in a quadratic free energy for -(/;_ . The free energy is minimized for 
tp- and the solution is used to get an effective free energy for ip^ . A subtle point is this derivation is the treatment of 
the relative phase between the two components ?/>+ and V_ . Within the approximations made below it can be shown 
that this phase is spatially uniform. The effective free energy for is then treated with standard methods to arrive 
at a London free energy. In particular, writing tp+ = \ip+ \ exp*^, taking \ip+\ to be spatially uniform and equal to the 
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homogeneous value, introducing V5 = V0 — ^A, keeping only terms of order v'^, minimizing this free energy with 
respect to A (to find as a function of b = V x A) yields the following London energy for vortex lattice in the A 
phase [the vortex lattice was introduced as a gaussian source (with a cutoff given by the A phase coherence length 
^A.+) for the London equation, as was done in Refs. P,p^] 



52 ^ e-'i'«^.+ 



^ — ^0 + TT^ T 2T2 — t" (3) 



where q is summed over all reciprocal lattice vectors of the vortex lattice, Bq is the spatially averaged magnetic 
induction, A^i = ^^^^ij'^ Ji-n^n penetration depth in the A phase when k — 0, and ^a- = "^^p^+^ft'K, ^® 

coherence length that describes the transition to the B phase from the A phase. In the B phase vortex lattice the 
London energy is 

q— ^-^ -D V K — K sin ^2; /t+At sin ^y/ B K'^ xi^ 

-yB hA^q2 J 

B,- 

where As = ^ aL-e^K penetration depth in the B phase when k — 0, simp — (note sintp — 1 for T = Tc^b), 

_ — ""^^^p^ is the coherence length that describes the transition to the A phase from the B phase, and £,b,+ is 
the coherence length associated with the B phase 5*+ order parameter component. We take kql — \b/^b.+ = 60 as 
is quoted by Huxley et al. |^ and for simplicity we take XbIS,b.+ ~ ^a/S,a.+ (which is correct to within 10%). 

The interesting physical behavior arises because ^ii'^'"^' = at the Tc b so that the non-local term (the term 

proportional to [k/k]^) becomes of order and thus as important as the local term. This is due to the divergence 
of the correlation length associated with the A to i? phase transition. Typically near Hd the only relevant length 
scale is the distance between vortices and the usual analysis can be applied However, near the A to i? phase 

transition the diverging correlation length appears as a second relevant length scale and, as will be seen below, plays 
an important role in the London theory. Corrections beyond the mean field model used here can be included in the 
London theory. This will change the temperature dependence of 1/Ca(b),-j but will not change the powers of q that 
appear in the non-local term [T^ ]. It is the q dependence of the London theory that is important in what follows. 
The complete London theory will also contain additional non-local terms (like those that appear in the borocarbides 
[10|), however these corrections are small because they formally carry an additional factor {£,a{B),+/^A{B))^ ~ 1/3600 
[10,||. Consequently, near Hd, the above London energy will describe UPts very well in the vicinity of the Tc.b- 
Below we discuss some observable consequences of this London free energy. 
Penetration Depth 

The temperature dependence of the non-local term becomes of order q^ at the Tc,b- This indicates that the 
temperature dependence of the penetration depth will exhibit a non-trivial behavior. For the field along the c-axis 
and a surface with a normal along the crystallographic a axis, the penetration depth in the B phase is given by 

= ; (5) 

l-fsin^+^ + ^(l-fsin^-^)^-4f^^ 

the result for the A phase is given by changing the label B to A and setting sinip — 1. This behavior differs from the 
usual Ginzburg Landau behavior. The penetration depth will also exhibit an in-plane anisotropy that is maximal at 
Tc_B. However, this anisotropy will be obscured by the existence of multiple domains of a symmetry breaking field. 
Vortex lattice structure 
The vortex lattice unit cell used here is defined in Fig. ^. 
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FIG. 2. Real space unit cell of the vortex lattice. 



The London energy derived above has the most significant impact on the vortex lattice structures for temperatures 
near Tc^b- EarUcr analysis implied a stretched hexagonal lattice near Tc^b which is not strongly field dependent 
(this was found for arbitrary k/k). Such a lattice results from a re-scaling of the variables qx and qy to remove a q^—qy 
orthorhombic symmetry term in the London energy | [Tl[ |. The resulting London energy in the re-scaled variables is 
isotropic (which implies a hexagonal vortex lattice in the re-scaled parameters). This analysis neglects the divergence 
of S,A{B).- the temperature approaches Tc^Bi which changes the results significantly. We performed a numerical 
analysis of the vortex lattice structure in the A phase for k/n = 0.36, the results are shown in Fig. ^. In this analysis, 
the q^ — qy anisotropy in the London energy is removed by re-scaling the q^ and qy pl] | . Since the re-scaled and original 
qx and qy differ by only a factor k/k, this re-scaling is ignored in the non-local term (note including the re-scaling 
in the non-local term is equivalent to keeping factors 0{k^ /k^) which have been previously thrown out). Fig. |^ only 
show results for H > 0.1\/2i?c (numerically it was found that Hd ~ 0.03^/2Bc note also Hc2 ~ 60\/2Bc)- As can be 
seen for T = Tc,b = 400 mK, the low-field vortex lattice phase is not hexagonal, but square for a substantial range 
of applied magnetic field. For T ^ Tc^b, the low-field vortex lattice is hexagonal. This can be understood by noting 
that £,A,- is finite for T ^ Tc.b and consequently, only the usual London theory applies at the lowest fields (since the 
vortices become infinitely separated as H Hd). However for T = Tc^b, S,a,- is infinite, so that the non-local terms 
in the London theory are important even as if — > Hd- Furthermore, it was found that for the lowest fields that were 
numerically accessible, the stable vortex phase for T = Tc^b was rectangular (a square lattice stretched along one of 
its sides). This lattice was found to be stable for Hd < H <sa 1.003-ffci (this field range is not shown in Fig. ||). In 
comparing the theory to existing experimental results, it is important to point out that the high-field vortex phase 
agrees with that observed by Huxley et al. |^ . We found that the high-field (nearly) hexagonal lattice has one of the 
basis vectors rotated 15 degrees from the crystallographic a axis which is in agreement Ref. |p|. This is a non-trivial 
result of the theory found here. In particular, the sign of the non-local term derived in Eq. ^ is opposite that used 
by Huxley et al. to explain the observed orientation. For smaller k/k, the vortex phase diagram is very similar to 
that found in Fig. |^ with one important difference: the region where the square vortex lattice phase appears becomes 
suppressed to lower fields and to temperatures closer to T^-^b- 
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FIG. 3. Vortex lattice structure as a function of applied field and temperature in the j4-phase for temperatures near the A 
to B transition {Tc,b = 440 mK). Results have been included for fields H > 0.1\/2Bc. The parameter C = 1 throughout the 
entire parameter range shown above. 



Given that the intriguing behavior found here was due the diverging correlation length at Tc.b, it is clear that 
this physics is not restricted to UPts. It will arise in other multi-component superconductors that exhibit multiple 
phase transitions. Some additional candidates for this behavior are Ui_a;Tha;Bei3 Sr2Ru04 p^-|l^], Celn3 [ p7| , 
electron doped cuprates and MgCNia [ p^pc| ]. In Ui_a;Tha;Bei3 multiple phase transitions have been observed 
[ p^ . Cubic Celna is likely to have a superconducting phase diagram similar to UPta with a two-component d-wave 
order parameter pl[ . Sr2Ru04 has been established as an unconventional superconductor with a two-component 
order parameter jljjl^ . The application of uniaxial stress should give rise to multiple transitions in this case . 
Electron doped cuprates have been argued to exhibit a d to s-wave pairing symmetry transition with doping ] . If 
this is indeed the case, then there should exist a particular doping where the two symmetries co-exist. This would 
imply multiple phase transitions. Finally, cubic MgCNia has been argued to be a spin-triplet superconductor for 
which it is likely that a multi-component order parameter arises [Q. In this case, the application of uniaxial pressure 
will then be required to give rise to multiple superconducting phase transitions. 

In conclusion, we have shown that for superconductors that exhibit multiple superconducting phases in zero applied 
magnetic field, the London theory in the vicinity of phase transitions takes on a very novel form. In particular, the 
diverging correlation length associated with the phase transition allows non-trivial vortex lattice structures to exist 
arbitrarily close to Hd- It was also found that the penetration depth will exhibit an unusual temperature dependence 
and anisotropy near such transitions. This theory was applied to UPts and the possible application of this theory to 
Ui_2,Th2;Bei3, Sr2Ru04, Celus, electron-doped cuprates, and MgCNia was discussed. 
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